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ABSTRACT 

This paper reports on the in-plane normal modes in the self-consistent and the cut-out 
power-law disks. Although the cut-out disks are remarkably stable to bisymmetric per- 
turbations, they are very susceptible to one-armed modes. For this harmonic, there is 
no inner Lindblad resonance, thus removing a powerful stabilising influence. A physical 
mechanism for the generation of the one-armed instabilities is put forward. Incoming 
trailing waves are reflected as leading waves at the inner cut-out, thus completing the 
feedback for the swing-amplifier. Growing three-armed and four-armed modes occur 
only at very low temperatures. However, neutral m = 3 and m — A modes are possible 
at higher temperatures for some disks. The rotation curve index (3 has a marked effect 
on stability. For all azimuthal wavenunibers, any unstable modes persist to higher tem- 
peratures and grow more vigorously if the rotation curve is rising (/? < 0) than if the 
rotation curve is falling {(3 > 0). If the central regions or outer parts of the disk are 
carved out more abruptly, any instabilities become more virulent. 

The self-consistent power-law disks possess a number of unusual stability properties. 
There is no natural time-scale in the self-consistent disk. If a mode is admitted at some 
pattern speed and growth rate, then it must be present at all pattern speeds and 
growth rates. Our analysis - although falling short of a complete proof - suggests that 
such a two-dimensional continuum of non-axisymmetric modes does not occur and that 
the self-consistent power-law disks admit no global non-axisymmetric normal modes 
whatsoever. Without reflecting boundaries or cut-outs, there is no resonant cavity and 
no possibility of unstable growing modes. The self-consistent power-law disks certainly 
admit equiangular spirals as neutral modes, together with a one-dimensional continuum 
of growing axisymmetric modes. 

Key words: celestial mechanics, stellar dynamics - galaxies: kinematics and dynamics 
- galaxies: structure - galaxies: spiral 



1 INTRODUCTION 



The preceding paper in Monthly Notices has presented the formidable machinery to permit direct solution for the 
in-plane normal modes of the family of razor-thin power-law disks. The results in this paper are based on a fast 
computer program that sets this machinery working. 

The entire spectrum of normal modes of a family of rigidly rotating stellar dynamical disks was deduced 



by Kalnajs (1972) in a bravura piece of mathematical analysis using "only pencil, paper and Legendre polyno- 
mials" (Toomre 1977). Differentially rotating disks are much harder and the results in the literature are sparse. 
Kalnajs (1977) devised a general procedure for deducing the normal modes of an axisymmetric stellar disk using ma- 
trix mechanics in action-angle coordinates and he subsequently applied his method to the isochrone disk ( pK^alnajs 1978 ) 
Over twenty years have passed since Zang's ( |1976D thesis on the scale-free disk with a completely flat rotation curve, 
yet it remains the most complete and instructive normal mode analysis available. In fact - perhaps daunted by the 
compl exities of mode analyses - many researchers hav e looked for convenient short-cuts, such as gas approxima- 
tions ( Bardeen 1978 ; Aoki et al. 197£ ; Bertin fc Lin 1996 ) or the use of cold disks with soften ed gravity ( Erickson 1974 ; 
Toomre 1981 ) or the combination of N-body work with approximate analytic treatments ( 3ellwood 1989 ). Our aim 
is to extend Zang's normal mode work to all the scale- free disks with rising and falling rotation curves, the power-law 
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disks. We are not alone in fixing attention on these disks as suitably simple models for which global stability analyses 
can be carried out. The gaseous power-law disks have attracted the attentions of a number of investigators (Schmitz 
& Eljert 1987|; |Lemos et al. 1991; Lynden-Bell & Lemos 1993; 3yer fc Tremaine 1996| ). The present paper gives 



a complete analysis of the global stability of the cut-out and the self-consistent power-law disks of stars. In the 
cut-out disks, stars close to the centre (and sometimes also at large radii) are prevented from participating in the 
disturbance. In the self-consistent disks, all the stars are mobile. 

Our study of the global stability of the cut-out power-law disks aims to determine how much random motion 
is required to stabilise a disk to modes of each angular harmonic. A recurring theme of numerical simulations 
throughout the seventies and eighties has been the ferocity of the bar instability , even in disks that are safel y 
stable to ax isymmetric perturbations ( Hohl 1971 ; Athanassoula fc Sellwood 1986 ; 3ellwood fc Wilkinson 1993 ) . 
Zang ( 1976 ), somewhat surprisingly, found that this is not the case for the scale-free disk with a flat rotation curve. 
It is bar-stable unless so cold that it is already susceptible to axisymmetric modes. For each power-law disk, it is 
interesting to locate the fastest growing modes and establish the physical origin of the instabilities. Are galaxies 
with rising rotation curves more or less stable than galaxies with falling curves? What kinds of density laws can be 
carved out to render a disk stable (or unstable)? As exact global stability analyses will remain scarce because of 
their cumbersome nature, it is important to assess how rules-of-thumb ~ such as Toomre's ( |1964| ) local criterion, or 
Ostriker and Peebles' ( 1973 ) global criterion - fare against exact results. 

The behaviour of the self-consistent disks must be very different from that of the cut-out disk. As originally 
pointed out by Kalnajs in 1974 (reported in Zang (1976)), all dependence on growth rate and pattern speed can be 



factored out of the integral equation for the normal modes. The same detail was noted independently by Lynden- 
Bell ^ Lemos (1993) using an argument based on dimensional analysis. Physical intuition suggests that as a disk 
is cooled below the critical temperature, growing modes become possible. But, the symmetry of the scale-free disk 
forbids a preference for one growth rate over another. At the critical temperature, the disk must admit modes with 
all frequencies or none. Thus, there seem to be two options: either there is a continuum of modes with all possible 
pattern speeds and growth rates, or there are no normal modes whatsoever. Do the scale-free, self-consistent power- 



law disks have normal modes? There are contradictory speculations on the answer to this question (Zang 1976 



Lynden-Bell fc Lemos 1993). 



The paper first concentrates on the cut-out disks and examines their global stability according to azimuthal 
wavenumber. Bisymmetric instabilities (m = 2) are studied in Section 3. The one-armed modes (m = 1) are 
examined in Section 4, while the other azimuthal harmonics (to = 3,4) in Section 5. Finally, Section 6 considers the 
self-consistent disk in the light of the insights gained in the study of cut-out disks. Before all this, however, we begin 
with a short introductory Section 2, which introduces some common notation and themes for the rest of the paper. 



2 THEORETICAL PREAMBLE 

Before discussing stability to any non-axisymmetric modes, it is helpful to start with a brief discussion of axisymmet- 
ric stability. This provides part of the intelle ctual framewor k needed to understand the non-axisymmetric results. 



For short wavelength, axisymmetric modes, Toomre (1964) showed that local stability requires that the velocity 



dispersion exceed a critical value du,iai-a- First, let us make the convenient abbreviation 

H5^- r[l(l-/3)]r[l(2 + /3)] 

"^^^^ r[i(i + /3)]r [1(2-/3)]' 

For the power-law disks, the critical velocity dispersion for stability to axisymmetric modes is 

3.3582830 

'^active fc\\ 

^" ^ = 27rL(/3)V2-/3 l^oq ' 

where Sactivc is the active surface density or the density spared by the cut-out function and E^q is the surface density 
implied by Poisson's equation. For the cut-out disks, au,min is a function of radius, although for the self-consistent 
disk it is not. We shall often relate the velocity dispersion du in both self-consistent and cut-out disks to the minimum 
velocity dispersion necessary to ensure local axisymmetric stability in the self-consistent disk. This ratio is referred 
to as Qs 



^ 3.3582830 ^ ' 

Local theory tells us that Qs > 1 is sufficient to stabilise against short wavelength axisymmetric modes. This result 



also holds good fo r stability against global axisymmetric modes, as has been found before by others ( Zang 1976 
Lemos et al. 1991 ) and is briefly demonstrated anew in Appendix |a]. 
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Table 1. This table shows the variation of some useful stability indicators with the rotation curve index (3 of the self-consistent power-law 
disks. Here, CTu min is the minimum velocity dispersion (measured relative to the local circular speed) needed to suppress axisymmetric 
Jeans instabilities entirely. X is the ratio of the circumferential wavelength of an m-armed disturbance to the longest Jeans unstable 
wavelength Acrit in a cold disk. 



Toomre (1964) showed that rotation alone stabihses axisymmetric modes with wavelengths in excess of Acrit 

_ AtT^GJ: _ 2nR Sactive , 

«2 - (2-/3)L(/3) S,, ■ 

Equivalently, lengthscales larger than Acrit are stable to Jeans instabilities even in a cold disk. It is often helpful to 
discuss non-axisymmetric instabilities in terms of a dimensionless ratio X (e.g., Toomre 1981). This is the ratio of 
the circumferential wavelength of an m-armed disturbance A to Acrit, i-e., 

^ A ^ 2tiR ^ {2~ (3)L{(3) Scg 

Acrit 'TlAcrit TTl Sactive 

Table 1 shows how these helpful stability parameters vary for the self-consistent power-law disks (Ecq = Sactivo)- 
The table makes the point that the properties of the disks change quickly as one moves away from the reference 
Toomre-Zang disk (/3 — 0). For example, Acrit increases by a factor of 1.75 on moving from (3 = io (3 — —0.5 and 
decreases by a factor of 1.64 on passing from /3 = to /3 = 0.5. The densities of any clumps formed simply through 
axisymmetric Jeans instabilities are almost a factor of three larger on moving from (3 = 0.5 to /3 = —0.5. Relative to 
the local circular speed, a-«,min roughly doubles or halves on making the same journeys. This table already suggests 
that the disks with /3 < will be harder to stabilise against Jeans instabilities and their kin than the disks with 
/3 > 0. The swing- amplifier is at its fiercest when \ < X < 2 and becomes weak once X > i (see Toomre 1981). So, 
Table 1 cautions us that the one-armed (m ~ 1) modes especially in the disks with rising rotation curves (/3 < 0) 
will be particularly troublesome to stabilise. 



3 BISYMMETRIC MODES 

This section provides the results of the stability analysis to global bisymmetric perturbations. From an observational 
standpoint, these are the most important, as some 30% of disk galaxies are strongly barred with perhaps a further 
30% showing evidence of weaker barring ( Bellwood fc Wilkinson 1993| ). There are even galaxies known which are 
unbarred in the optical but are barred in the infrared. Approximate bisymmetry is a characteristic of most disk 
galaxies. 



3.1 The Behaviour of the Eigenvalues 

3.1.1 The Effect of Growth Rate and Pattern Speed 

Let us recall that the mathematical eigenvalue A is the ratio of the response to the imposed density transforms 
(see Section 4.3 of the preceding paper). Mathematically, it is defined by Ares (a) = AAimp(a). If the modulus of 
the complex mathematical eigenvalue |A| is greater than unity, the response is more vigorous than the imposed 
perturbation. If |A| < 1, the response is less than the original disturbance. When the mathematical eigenvalue is 
unity, the response is identical in magnitude and phase to the imposed disturbance. This is a self-consistent mode. 
To discover global modes, it is important to understand the position of the eigenvalue in the complex plane and to 
learn what choices of growth rate, pattern speed and temperature can bring it to the point (1,0). Cool disks are 
(usually) less stable than warm ones. By increasing the temperature sufficiently, the modulus of the mathematical 
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13= 0.?5 ^ = -0.?5 

Figure 1. The dependence of the largest mathematical eigenvalue on growth rate and pattern speed for bisymmetric disturbances in 
a irmer cut-out disk with N = 2, Qs = 1 and /3 = ±0.25. [The solid lines show curves of constant s at intervals of 0.1 from s = 0+ 
to s = 0.5. The dotted lines show curves of constant Qp at intervals of 0.1 from Qp = 0.5 to Hp = 2.0. For fS = 0.25, ct^ = 0.283 and 
7 = 11.0; for /3 = -0.25, an = 0.509 and 7 = 3.36]. 



eigenvalue shrinks so that no choice of growth rate s and pattern speed Sip enables it to attain the value (1,0). On 
physical grounds, slowly-growing disturbances are expected to be more easily excited than faster ones. Thus, the 
modulus of the mathematical eigenvalue is greater for smaller growth rates. The pattern speed has a twofold effect, 
as it controls the phase of the response and determines the position of the co-rotation radius Rqr and the inner and 
outer Lindblad resonances, i?iLR and i?oLR- For the power-law disks, these crucial radii occur at: 



RiLR I m- \ i?OLR _ / m + V2~^ y+" i^cR ' ^ ^ 

Rq \ mQp I Ro \ mClp I ' -Ro 




(6) 



These formulae are given for general aziniuthal wavenumber m. Henceforth, the emphasis in this section is on m = 2. 
Any modes are generally confined to the region between the two Lindblad resonances. For Op > 1 + 1^2-/3, the 
outer Lindblad resonance lies inside the inner cut-out a,t R = Rq. The disturbance is restricted to a part of the disk 
where there is not much matter free to respond. The re sponse , and hence the mathematical eigenvalue, is expected 
to be small for large pattern speeds. For flp < 1 — i\/2 — /3, the inner Lindblad resonance lies beyond the inner 
cut-out. Now, the presence of the inner cut-out has a minimal effect on the response, since the density experienced 
is similar to that in the self-consistent disk. Each imposed logarithmic spiral then excites only a pure logarithmic 
spiral in response and so the eigenvalue is expected to be real. In a doubly cut-out disk, with sufficiently low pattern 
speed, the inner Lindblad resonance can fall beyond the outer cut-off at i? = i?c. In this case, the mathematical 
eigenvalue is again expected to be very small. 

Fig. H shows the effect of growth rate and pattern speed on the mathematical eigenvalue for = 2 disks with 
/3 = ±0.25. The temperature of each disk is chosen so that it is locally stable to axisymmetric perturbations, i.e., 
Qs — I- The curves have the general behaviour expected from our theoretical discussion. As the pattern speed is 
increased beyond ~ 1, the modulus of the mathematical eigenvalue declines. For large f2p, the eigenvalue turns 
towards the origin and the growth rate becomes smaller. The disk responds to the perturbation rather weakly, since 
most of the mass has been removed from the central region. Now, the mathematical eigenvalue depends on growth 
rate s and pattern speed f2p through the single complex frequency lo = mfip -I- is. Since the kernel is analytic, the 
mathematical eigenvalue A is an analytic function of uj. As realised by Zang ( 1976| ), this means that X{uj) gives 



a conformal mapping of the upper half of the complex cj-plane (since the growth rate is positive) to the complex 
A-plane. Conformal mappings preserve the angles of infinitesimal polygons. So, curves of constant s and constant 
flp are orthogonal. Fig. |l| shows a problem which can occur. For /3 = 0.25, there is a discontinuity near (fip — 0.5, 
s = 0.4), even though A is expected to change smoothly. The discontinuity arises because here the second largest 
eigenvalue overtakes the previously largest one, so that the same eigenvalue is not being followed. 

What are the values of the growth rate, pattern speed and temperature for which the eigenvalue is unity? Fig. |^ 
already suggests that, for Qs — I disks with /3 in the range —0.25 to 0.25, no choice of growth rate and pattern speed 
results in an eigenvalue equal to unity. For each /3, the curves of constant growth rate cross the real axis well to the 
left of unity. From our theoretical discussion, the mathematical eigenvalue is likely to continue towards the origin as 
the pattern speed increases still further. Assuming that, at any pattern speed, the marginal eigenvalue curve has the 
largest modulus and that the curves of constant s do not cross, it seems improbable that any solutions can exist at 
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Figure 2. Marginal eigenvalue curves for bisymmetric disturbances in inner cut-out disks with N = 2 and values of the rotation curve 
index from /3 = —0.5 to /3 = 0.35 in steps of 0.05. The velocity dispersion for each /3 is that corresponding to local axisymmetric stability 
{Qa = 1). [The plot shows the largest mathematical eigenvalues for vanishing growth rate, s = lO"'^, and 20 values of the pattern speed: 
ftp = 0.0001, 0.0005, 0.001, 0.005, 0.01, 0.02, 0.05, 0.1, 0.2, 0.3, . . . 1.2, 1.3.] 

pattern speeds lower than those shown. Fig. |l| plots the behaviour of the dominant eigenvalue, i.e., that which has 
the largest modulus at a particular pattern speed and growth rate. For the range of pattern speeds investigated so 
far, the modulus of the dominant eigenvalue has not reached unity, and so none of the lower eigenvalues can yield 
a mode. Let us also observe from Fig. |] that as /? is reduced from 0.25 to —0.25, the set of curves tilts upwards, 
crossing the real axis at higher values, and thus coming closer to possessing a unit eigenvalue. The disks with rising 
rotation curves (negative (3) are less securely stable than their relatives with falling rotation curves (positive f3). 
Fig. H shows marginal eigenvalue curves for several different values of (3 in the range —0.5 to +0.35. The curves 
are drawn for Qs = 1, so all the disks shown are already stable to axisymmetric perturbations. The disks with 
rising rotation curves have higher velocity dispersions than those with falling rotation curves, but even so they are 
more prone to bisymmetric disturbances. As /3 is decreased, the magnitude of the marginal eigenvalue at a given 
pattern speed increases. Eigenvalues in disks with rising rotation curves cross the real axis at higher values than for 
their cousins with falling rotation curves. For values oi f3 < —0.470, Qs = 1 is insufficient to stabilise the disk to 
TO = 2 perturbations. Disks whose rotation curves rise faster than this succumb to bar-like modes, even if stable to 
axisymmetric modes. 

3.1.2 The Effect of the Temperature and the Inner Cut- Out 

Raising the velocity dispersion (t„ moves the marginal stability curves closer to the origin, while lowering it moves 
them further out. At some value a^u^min, the marginal stability curve will pass through (1,0); this is the temperature 
at which the disk is just stable to bisymmetric perturbations. Fig. || shows how the temperature affects the marginal 
stability curve. The eigenvalue depends separately on the pattern speed and temperature, so the mapping from the 
(rjp-Qs) plane to the complex A plane is not conformal. This is evident from Fig. ^, where lines of constant flp and 
constant Qs do not intersect at right-angles. Again, this demonstrates that power-law disks with falling rotation 
curves are securely stable to bisymmetric perturbations. Even when the velocity dispersion is only 60% of that 
needed for axisymmetric stability, the (3 — 0.25 disk admits no to = 2 modes. 

The inner cut-out strongly influences the stability of the disk (c.f., Zang 1976; Toomre 1977). To illustrate this. 
Fig. ^ shows the marginal stability curves in disks with different /? and different values of the inner cut-out index N. 
Increasing the cut-out index N makes the disk much more susceptible to instabilities. Curves with higher N cross 
the real axis at successively higher values of A. Disks with iV = 1 are highly stable to to = 2 perturbations. To 
understand why the inner cut-out has such a profound effect on the stability of the disks, let us consider the position 
of the Lindblad resonances. For low values of the pattern speed, the inner Lindblad radius lies outside the inner 
cut-out radius Rq. Incoming waves are absorbed by the inner Lindblad resonance, which thus damps disturbances 
in the disk and reduces the possibility of a mode. However, once fip exceeds ~ 0.25, the inner Lindblad resonance 
moves inside the inner cut-out radius. If the cut-out is gentle {N =1), the incoming waves pass through it, reach the 
inner Lindblad resonance and are absorbed. If however the cut-out is sufficiently sharp (as for larger cut-out indices), 
it presents a barrier which reflects the incoming trailing waves. Instead of reaching the inner Lindblad resonance and 
being absorbed, they are reflected back as outgoing leading waves. The exceptional stability of iV = 1 disks is mostly 
due to the fact that the surface density is disturbed much less abruptly than for higher cut-out indices, presenting 
less of a reflective barrier to incoming waves. Additionally, there is a much greater surface density within the inner 
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Figure 3. The dependence of the largest mathematical eigenvalue on temperature and pattern speed for inner cut-out disks with 
/3 = ±0.25. Curves are plotted for bisymmetric perturbations with vanishing growth rate in a disk with cut-out index N = 2. The 
temperature used for each curve is expressed in terms of Qs. [For fS = 0.25, 0"u = 0.283Qs; for /3 = —0.25, CTu = 0.509Qs]. 





/3 = 0.25 (3 = -0.25 

Figure 4. The dependence of the largest mathematical eigenvalue on pattern speed, for different inner cut-out indices A'^. Curves are 
plotted for bisymmetric perturbations with vanishing growth rate in an inner cut-out disk with /3 = ±0.25. The curves are labelled with 
the value of the inner cut-out index A'^. [In each case Qb = 1. The diamonds show the results for f2p = 0.1, 0.2, 0.3 ... 2.O.] 



cut-out when iV = 1; the active surface density diverges at the origin when the rotation curve is faUing, and even 
when the rotation curve is rising, the surface density remains large to much smaller radii. By contrast, for > 1, 
the active surface density rapidly falls to zero within Rq. Waves are thus suppressed by simply having no medium 
through which to propagate. 



3.2 The Marginal Modes 

The marginal modes are those with vanishingly small growth rate. How cold does the disk have to be for this 
transition from stability to instability to occur? This question can be answered by finding the temperature for which 
the marginal eigenvalue curve passes through (1,0). Fig. ^ shows the minimum temperature needed for bisymmetric 
stability as a function of /3. The solid lines show the results for the inner cut-out disks with iV = 2, 3, and 4. The 
broken lines show results for the corresponding doubly cut-out disks with AI = 2, 4, 6 and — 10. Curves with 
different N lie well apart. The more sharply the core is cut out, the more unstable the disk becomes. For each N, 
the curves with different M lie very close together. Even with this low truncation radius, the outer cut-out function 
makes little difference to the temperature necessary for stability. With a larger Rc, the doubly cut-out curves are 
scarcely distinguishable from those with only an inner cut-out. The right-hand plot in Fig. ^ can be used to compare 
the relative stability of the various cut-out disks to m = and m = 2 perturbations. The N = 2 disk is less prone 
to bisymmetric than to axisymmetric instabilities over almost all f3. For iV = 3, the relative tendency to m = 2 and 
m = disturbances depends on (3. Roughly speaking, iV = 3 disks with rising rotation curves are more susceptible 
to m = 2 disturbances whereas disks with falling rotation curves are more susceptible to to = 0. Conversely, disks 
with iV = 4 are much more prone to bisymmetric disturbances over almost all /?. Only when the rotation curve is 
steeply falling does the to = mode become harder to stabilise. It is perhaps surprising that it is the lowest outer 
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Figure 5. Minimum temperature for bisymmetric stability plotted against rotation curve index /3 for disks with various cut-out functions. 
The left-hand plot shows the minimum velocity dispersion cru\ the right-hand one the same data presented in terms of the stability 
parameter Qs- The solid lines indicate the results for inner cut-out disks with A'^ = 1, 2, 3, 4. The broken lines indicate the corresponding 
doubly cut-out disks with M = 2, 4, 6 and Rc = 10. 



^ N=2 




-0.4 -0.2 0.0 0.2 0,4 

Figure 6. Critical pattern speed Hp plotted against rotation curve index j3 for m = 2 modes in disks with various cut-out functions. 
The solid lines are labelled with the inner cut-out index N. For each A^, different outer cut-out indices M are plotted. 



cut-out index (Af — 2) which shows the largest departure from the inner cut-out curves. In this instance, it seems 
that tapering the disk gently at large radius has more of an effect than truncating it abruptly. This is presumably 
because for lower M, the effect of the outer cut-out is felt further in towards the centre of the disk. For higher M, the 
outer cut-out is sharp, but occurs entirely outside the outer Lindblad resonance. Disks with a gentle outer cut-out 
are more stable than those where the cut-out is sharper. The critical pattern speed obtained for these marginal 
modes is shown in Fig. ^ For all cut-out functions and /3, the pattern speed is confined to a fairly narrow range, 
roughly 0.3 to 0.6. Over a wide range in /3, the pattern speed varies roughly linearly with (3. 

Let us look at the shapes of some of these marginally stable modes. The contour plots in Fig. ^ show the form 
of the modes for (i = ±0.25 in inner cut-out disks with N = 2. For falling rotation curves, the modes are more 
tightly wound than for rising rotation curves. They extend beyond co-rotation up to the outer Lindblad resonance, 
whereas in disks with rising rotation curves the modes are concentrated between the inner Lindblad resonance and 
the CO- rotation radius. The modes are also more tightly wound for disks with gentle cut-outs. These effects are 
caused by the higher temperatures needed for marginal modes when the rotation curve is rising. If a growing mode 
exists, then the stars do not have enough random motion to move entirely out of a density enhancement before it 
has grown significantly. As the temperature increases, the stars have more random motion, and the spiral modes can 
therefore be more loosely wound. As shown in Fig. the pattern speed of modes increases with f3 and is greater 
for disks with lower cut-out indices. Since the Lindblad resonances move inwards as the pattern speed increases, 
the pattern is concentrated within a smaller radius for /3 — 0.25 than for (3 = —0.25, and similarly for the lower 
cut-out indices. Fig. ^ shows the position of the Lindblad and co-rotation resonances for modes in cut-out disks. 
It is significant that, for every value of f3, the pattern speed is large enough to push the inner Lindblad resonance 
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Figure 7. Density contour plot of a marginally stable mode in inner cut-out disks with N = 2. The solid lines mark 10%, 20%, 40%, 60% 
and 80% of the maximum density. The dotted lines show the nodes, where the perturbation density vanishes. The Lindblad resonances 
and the co-rotation radius are indicated by dashed circles. [The left-hand plot is for /3 = 0.25, for which the marginally stable mode has 
Qs = 0.500, 7 = 48.4, ct„ = 0.142, Hp = 0.596 and Qu 20.4. The right-hand plot is for /3 = -0.25, for which the marginally stable 
mode has Qs = 0.631, 7 = 9.19, a„ = 0.321, Qp = 0.533 and 8.4]. 
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Figure 8. Position of resonances for marginal m = 2 modes in cut-out disks with different rotation curve index /3. The dashed horizontal 
lines show the position of the inner and outer cut-out radii. The pattern speed corresponding to the marginal mode in each disk (Fig. ^ 
is used to calculate the the resonance radii. The solid lines mark the Lindblad resonances; the dotted line marks co-rotation. [The three 
plots show results for N = 2, 3 and 4; each plot shows data for doubly cut-out disks with M = 2, 4, 6, Rc = 10]. 



just inside tlie inner cut-out. Tlie disk does not admit slowly-rotating growing modes for which the resonance would 
emerge beyond the cut-out. The resonance is allowed nearer the cut-out when the cut-out is sharper, and hence a 
more efficient reflector. Similarly, on examining modes in an unstable disk (Read 1997), the fastest-growing modes 
have the largest pattern speeds. As the pattern speed decreases, allowing the inner Lindblad resonance to move out 
closer to the cut-out radius, the growth rate slows. 



3.3 The Ostriker-Peebles Criterion 



Ostriker and Peebles ( 1973 ) suggested that disk galaxies are stable to bar- like modes only when the ratio of total 
rotational energy T to total gravitational energy \W\ is less than 0.14. The energies W and T for the self-consistent 



Stability of Power-Law Disks II - The Global Spiral Modes 9 




Figure 9. Comparison of the results for the cut-out disks with the Ostriker-Pcebles criterion. The plot shows the ratio of total rotational 
energy to total potential energy necessary for marginal stability to bar modes. The solid lines are labelled with the appropriate inner 
cut-out index N. For each N, different outer cut-out indices M are plotted, as shown in the key. 



disk are given in Section 2 of Paper I. We thus deduce the ratios as ( Evans 1994 ) 
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result can be derived from the values for /3 7^ using a property of the gamma function ( prad- 



Note that the /3 

shtey fi fc R}^k 1978| , 8.328.2). 

Fig. |9| shows the ratio T/|T/F| calculated for the self-consistent disk (eqs. (7) and (8)). The anisotropy parameter 
used at each /3 is that needed for marginal stability to bisymmetric modes in the cut-out disk. Note that can 
exceed 0.5, since the standard form of the virial theorem does not apply to the power-law disks, which are of infinite 
extent. For all the disks, the ratio is much greater than the value of 0.14 suggested by Ostriker and Peebles. 

This demonstrates that the Ostriker-Peebles criterion is not useful in predicting the stability of the power-law disks. 
However, the temperatures used in Fig. ^ are derived for the cut-out disks. As the self-consistent and cut-out disks 
have different stability properties, perhaps the value of should be calculated for the cut-out disks. In this 

case, the density S and streaming velocity (v) must be computed numerically using the cut-out distribution function. 
Details of this calculation are given in Read (1997) - there is only a small change in the ratio T/|iy| and the main 
conclusion of this Section is unaffected. For the power-law disks, the Ostriker-Peebles criterion is untrustworthy as 
it greatly overestimates the amount of energy in the form of random motion necessary to achieve stability against 
bar modes. 



4 ONE-ARMED MODES 



The observational motivation for studying one-armed modes is that mild lop-sidedness is seen in many spiral galaxies, 



with substantial asymmetries present in some late- type galaxies (Baldwin et al. 1980). This is evident in both the 
optical and HI data (Richter & Sancisi 1994; Rix & Zaritsky 1995). One-armed disturbances differ from many-armed 



disturbances in one crucial respect. It is intuitively clear that one-armed patterns can move the barycentre from 
the origin towards the overdense spiral arm. For an isolated system, acceleration of the barycentre is not physical. 
The possibility occurs here because the analysis holds the equilibrium disk fixed while imposing a perturbation 
on it. If one-armed waves grow in a galaxy, the equilibrium disk is displaced to balance the perturbation so that 
the ba rycentre rema ined fixed. Our analysis thus reproduces a famous mistake by Maxwell in his Adams Prize 
essay ( Maxwell 1859 ). A complete study of the m — I modes requires that our analysis be modified to permit the 
barycentre to move. Instead, we investigate the stability of the disk to "nai've" one-armed disturbances and check a 
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Figure 10. The dependence of the largest mathematical eigenvalue on growth rate and pattern speed for one-armed disturbances in 
inner cut-out disks with N = 2 and Qs = 1. The horizontal axis is Re[A] and the vertical axis is Im[A]. [In the plot, /3 = 0.25 (cti, = 0.283 
and 7 = 11.0). The solid lines show curves of constant s at 11 values oi s: s = 10"'^, 0.0001, 0.0005, 0.001, 0.005, 0.01, 0.05, 0.1, 0.2, 0.3, 
0.4. The dotted lines show curves of constant Hp at 20 values of Q.^: = 0.0001, 0.0005, 0.001, 0.005, 0.01, 0.02, 0.05, 0.1, 0.2, 0.3 .. . 
1.2, 1.3]. 



posteriori that the barycentre moves only shghtly. In this, we are in the good company of Toomre and Zang, who 
did exactly this in their analysis of the disk with a flat rotation curve. 



4.1 The Behaviour of the Eigenvalues 

Let us first investigate how the largest mathematical eigenvalue depends on pattern speed and growth rate for the 
temperature at which the self-consistent disk is locally just stable to axisymmetric disturbances. Fig. |l^ shows the 
eigenvalue curves for a disk with j3 = 0.25 and N = 2. For growing perturbations (s > 0.05), these curves are 
qualitatively similar to those obtained for bisymmetric perturbations (Fig. |l|). Curves of a given s cross the real 
axis at higher values for these m = 1 plots than for the earlier m = 2 plots, indicating that the disk is more prone 
to one-armed than two-armed disturbances. For instance, we saw in Fig. |^ that Qs — ^ was already sufficient to 
stabilise the f3 — 0.25 disk to m = 2 disturbances. But, Fig. ^ indicates that this disk remains susceptible to m = 1 
disturbances. One such mode is indicated by the dashed line. The most striking difference between the present 
situation and that for m = 2 occurs as the growth rate is reduced to zero. It appears that in the limit of vanishing 
growth rate, the marginal eigenvalue curves never cross the real axis. Eigenvalues for small but non-zero growth 
rate (e.g. s — 10~^, 10""^) follow the marginal curve for high pattern speeds. As the pattern speed is reduced, one 
by one they peel away from the marginal curve and cross the real axis. This has dire consquences for the stability 
of the disks to one-armed disturbances. Ordinarily, disks may be stabilised by raising their temperature, but this 
does not work here. As the temperature is increased, the curves with higher growth rate do indeed shrink towards 
the origin. But the curves with very low growth rate continue to stretch out along the real axis. It appears that, no 
matter how high the temperature, a self-consistent mode exists for sufficiently low growth rate and pattern speed. 
Unlike for m — 2, there is no clear critical temperature distinguishing stable and unstable disks. We can only find a 
temperature sufficient to stabilise the disk down to \u!\ — 10^"^, lO^"' or whatever. Thus, it seems that there is always 
at least one mode. The existence of further modes depends on the behaviour of the sub-dominant eigenvalues. The 
sub-dominant eigenvalues behave in a similar way to the dominant eigenvalue as the growth rate and pattern speed 
are varied. 

What is the physical mechanism underlying the instability? A distinguishing feature of one-armed perturbations 
is that they are the only non-axisymmetric disturbances which have no damping inner Lindblad resonance. In the 
absence of an inner Lindblad res onance, incoming trailing waves are reflected b y the inner cut-out as leading waves, 



which are then swing-amplified ( Goldreich fc Lyndcn-Bell 1965 ; Toomre 1981 ). For m ~ 2, many of the incoming 



waves are absorbed by the Lindblad resonance and never return from the central regions. This point is illustrated in 
Fig. These plots show the density transforms of m = 1 and m = 2 modes in disks of various temperatures, chosen 
so as to obtain modes of broadly similar growth rates. The top plot in each figure shows a quickly-growing mode 
in a cool disk. Both transforms arc dominated by a large trailing component. Heating the disk reveals differences 
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Figure 11. Density transforms for m = 1 (left, = 0,25) and m = 2 (right, /? = —0.25) modes in disks of various temperatures. The 
captions to each plot record the value of Qs, and the growth rate and pattern speed of the mode. For m = 2, the position of the inner 
Lindblad resonance is also stated. In each plot, the solid curve shows Re[j4(a)], the dotted line Im[A(a)], and the dashed line the envelope 
±|j4(a)|. The dashed vertical line separates leading from trailing waves. 



between the m = 1 and m = 2 modes. For m = 1, the presence of a large leading component - previously hidden 
by swing-amplified trailing waves - becomes apparent. Close to marginal stability, in the bottom plot, the leading 
and trailing components are of almost equal amplitude. However, for m = 2 modes, no such leading component 
exists. This indicates that few leading waves are returning from the central region to fuel the amplifier. In fact, the 
physical origin of the m = 2 modes is quite diflferent. They occur in such cold disks that they are basically just Jeans 
instabilities. 



12 A^. W. Evans and J. C. A. Read 




Figure 12. Growth rates of the four fastest-growing one-armed modes for inner cut-out disks with Qs = 1-0. The four plots show results 
for Af = 1, 2, 3, 4. 




Figure 13. Pattern speeds of the four fastest-growing one-armed modes for inner cut-out disks with Qs = 1.0. The four plots show 
results for Af = 1, 2, 3, 4. 



4.2 Growing Modes 

How does the rotation curve and inner cut-out affect the pattern speed and the growth rate of the one-armed 
disturbances? To answer this, we need a way of comparing the stabihty of different disks. In the m = 1 case, there is 
no critical temperature above which the disks are stable. So, instead we employ the growth rate of the fastest-growing 
mode in disks with the same temperature as a measure of stability. More unstable disks admit faster-growing modes. 



Fig. 12 shows the growth rates obtained for the four fastest-growing modes in four different inner cut-out disks 
with Qs = 1- Fig- |lj shows the corresponding pattern speeds. Even though the condition Qs = 1 corresponds to 
larger velocity dispersions for negative (3, modes still grow faster in disks with rising rotation curves. Disks with 
sharper cut-outs are more unstable than those where the cut-out is gentle. For one-armed modes, there is no inner 
Lindblad resonance, so even disks with iV = 1 admit growing modes. However, it appears that a sharper cut-out still 
destabilises the disk, presumably because the increased reflection makes a more efficient feedback circuit for swing 
amplification. 

Figs. show contour plots of the fastest-growing modes in inner cut-out disks with Qs~^- The dominant 

wavenumber of the modes - that is, the wavenumber at which the magnitude of the density transform is greatest 
- is reported in the caption of each figure. As for the marginally stable bisymmetric modes plotted in Fig. 0, the 
one-armed spirals shown here are more tightly wound in disks with falling rotation curves than disks with rising 
rotation curves. In accord with the lack of an inner Lindblad radius, these one-armed modes extend right to the 
centre of the disk. Modes in disks with rising rotation curves have more density concentrated close to the centre than 
those in disks with falling rotation curves. When discussing the distance to which these modes extend out in the disk, 
it is useful to distinguish between the distance in physical space, and the distance in terms of the co-rotation radius. 
Modes in disks with low velocity dispersions generally extend further beyond the co-rotation radius than modes in 
warmer disks. However, modes in cooler disks also rotate more quickly, so their co-rotation and outer Lindblad radii 
are smaller. Thus the actual extent in space of a mode is generally smaller in cool disks than in warm disks. Bearing 
these temperature effects in mind, there is little difference in terms of spatial extent between one- and two-armed 
modes. The one-armed modes shown in Figs. ^-^| extend further in space, but less far beyond the co-rotation 
radius, than the two-armed modes plotted in Figs. |7|. This reflects the higher temperature of the disks in which the 
m = 1 modes are plotted. 



4.3 The Position of the Barycentre 

One-armed modes are unique in that they are capable of displacing the barycentre of the disk away from the origin. 
Visual inspection of Figs. |l^ and ^ suggests that in these modes the barycentre is still close to the origin. The spiral 
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/S= 0.25 |S = -0.25 

Figure 14. Density contour plot of the fastest-growing m = 1 mode in inner eut-out disks with A'^ = 1 and Qs = 1. The soUd lines mark 
10%, 20%, 40%, 60% and 80% of the maximum density in this range; the dotted lines show the nodes, where the density perturbation 
is zero. The dashed eircle indicates the co-rotation radius. [The left-hand plot is for (3 = 0.25, for which the fastest-growing mode has 
f2p = 0.0974, s = 0.0140, Ou Ri 4.2. The right-hand plot is for /3 = -0.25, for which the fastest-growing mode has Hp = 0.0807, s = 0.0343, 
au » 2.8]. 




Figure 15. Density contour plot of the fastest-growing m = 1 mode in inner cut-out disks with N = 4 and Qs = 1. [The left-hand plot 
is for /3 = 0.25, for which the fastest-growing mode has Op = 0.214, s = 0.100, au ~ 4.6. The right-hand plot is for /3 = —0.25, for which 
the fastest-growing mode has Clp = 0.177, s = 0.139, au RJ 2.8]. 
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Figure 16. The displacement of the barycentre from the origin for A'^ = 1 disks with (3 = ±0.25. The displacement is plotted in units of 
the co-rotation radius. We have set Ep = So to make the point that the barycentre does not move to second order. The dashed vertical 
lines mark the position of the co-rotation and outer Lindblad radii. 



arm wraps around the origin, with more density in the parts of the arm closer to the origin. It is difficult to guess by 
eye even the direction of any displacement of the barycentre from the origin. Intriguingly, the density distribution 
in one of the cases shown has two separate peaks, one on either side of the origin. This further suggests a balancing 
of the density on either side of the origin. 

In general, the coordinates of the barycentre of the part of the disk within a radius L are given by 



M{L)x 



dR / deT.{R,e)R^ cosO, 



M{L)y 



dR / de'E{R,0)R'^ sine, 



(10) 



where M{L) is the mass contained within the radius L. Each logarithmic spiral component makes equal positive 
and negative contributions at every radius, so the mass M{L) contained within a radius L is just the integral of the 
equilibrium density. Conversely, as the equilibrium density is axisymmetric, it does not shift the barycentre from the 
origin. So, the expression (10) for the x-coordinate of the barycentre becomes (using eqs. (44-45) of Paper I) 
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where we must take the real part to obtain the physical co-ordinates. For m ~ 1, we obtain 



l-/3Ep 



2 So \i?o 



Re 



+ 00 



daA{a) 



Rr. 



(12) 



--fine 



2 So 



L 
i?o 



1/2+13 



Im ■ 



+ 00 



daA(a) 



L 
Ro 



(13) 



If the integral over a is non-zero, the barycentre is displaced from the origin. As L is taken to infinity, the coordinates 
of the barycentre tend to infinity for /? > — ^. If we consider a large but finite value of L, eqs. ^) and (|l|) show 
that the barycentre gradually spirals outwards. The distance of the barycentre from the origin is given by 
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Fig. |l^ shows the distance from the origin of the barycentre of the part of the disk within a radius L in units of 
the co-rotation radius, i.e., \J -\- / Rck- The co-rotation radius roughly delimits the edge of the spiral mode, so 
this quantity describes how far the barycentre is displaced from the origin relative to the overall size of the mode. 
The dashed vertical lines in Fig. |l^ indicate the positions of the co-rotation and outer Lindblad resonances. For 
small values of L, the displacement of the barycentre is approximately linear. This is as expected - if we consider 
only the very inner region of the disk, the emerging spiral arm is predominantly on one side of the origin. But as 
we move outwards, the contributions of the spiral arm from opposite sides of the origin begin to cancel out. For 
values of L beyond the outer Lindblad resonance, the part of the disk under consideration is now large enough to 
include the entire spiral mode. The forces pulling the barycentre away from the origin in different directions cancel 
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almost exactly, so that the barycentre of the whole disk remains close to the origin. The fast oscillation of the 
term exp(ialnL) is probably the reason for the cancellation. For typical modes, A{a) is a smooth function of a 
which tends to zero as a —^ ±00. As L is increased, the term exp(zalnL) must eventually vary much faster with 
a than does A{a) (unless A{a) is a delta-function or some other strongly peaked entity). Thus the integration over 
wavenumber in eq. ( p^ adds up many positive and negative contributions of nearly equal magnitude. These roughly 
cancel, ensuring that the integral almost vanishes. At least in the linear regime, it appears that self-consistent modes 
do not shift the barycentre of the disk substantially. 



5 THE HIGHER ANGULAR HARMONICS 



This section briefly examines modes with m = 3 and 4. These are dubbed respectively triskele and tetraskele modes, 
after the three- and four-legged devices common in early European art. The observational motivation for studying 
these modes comes from Fourier decompositions of the spiral patterns of real galaxies, which indicate significant 
contributions from m = 3 and m — 4 harmonics (Elmegreen & Elmegreen 1985| ). The same result is suggested by 
analysis of N-body bars, which tend to be more rectangular than purely elliptical (Sparke & Sellwood 1987). 



5.1 A Numerical Delicacy 

Modes with m > 2 present an additional numerical delicacy. The expressions for the angular momentum function 
derived in Appendix C of Paper I involve division by the combination Ik + mil. When Ik + mn vanishes, the integral 
still exists in a Cauchy principal value sense. In the modes examined so far, this problem has not arisen. For m = 0, 
lii + mD, vanishes at the I — radial harmonic; but in this case the angular momentum function Fqq is simply zero. 
For m — 1 and m = 2, Ik + mfi never vanishes. But for higher angular harmonics, such as m = 3 and m ~ 4, the 
term Ik + mfl may vanish at some eccentric velocity. 

What does it mean physically when Ik + 111^1 vanishes? The time the star takes to perform / complete revolutions, 
27r|Z|/r2, is then identical to the time it takes to perform m radial oscillations, 27rm/K. Thus the star's orbit closes 
in an inertial frame. Closed orbits such as these do not occur in all disks. For Ik -\- mfl to vanish, then in terms of 
the auxiliary integral (defined in eq. (18) of Paper I), J2{U) = — 27rZ/m. In section 2 of Paper I, we saw that J2{U) 
declines monotonically from 27r/-y/2 — (5 a,t U = to tt at [/ = 00. Thus Ik + m(l will vanish for some finite U if 
there is an integer / which satisfies 

Whether or not this occurs depends on the values of m and /?. For m — 3, Ik + mCl vanishes at I — —2 for (3 > —0.25, 
while for /3 < —0.25 it never vanishes. For (3 = —0.25, Ik + mCt vanishes at U — 0, meaning that the relevant orbits 
are circular. For m = A, Ik mQ vanishes at I = —3 for /3 > 2/9, while for (3 < 2/9 it never vanishes. Even when 



Ik + mfl vanishes, the mathematical eigenvalue can still be evaluated accurately (Read 1997). The Gauss-Laguerre 
quadrature used for the integration over eccentric velocity successfully interpolates across the troublesome value of 
U. 



5.2 Triskele or Three- Armed Modes 

5.2.1 Global Stability of the Cut- Out Disks 

Fig. |l^ shows the dependence of the largest eigenvalue on growth rate and pattern speed for /? — 0.25. The left 
panel shows the results for Qs = 1.0, the right panel for Qs = 0.5. This plot already suggests that the power-law 
disks are very stable to m = 3 perturbations. Even when the temperature is as low as Qs = 0.5, the eigenvalue 
curves intersect the real axis at values considerably less than unity. For disks with N = 2, Qs has to be lower than 
0.4 before even the (3 — —0.5 disk becomes unstable to m = 3 perturbations. For disks with falling rotation curves, 
even lower temperatures are necessary. At these very low temperatures, the convergence of the numerical method is 
somewhat untrustworthy. This is mainly because the modes are more tightly wound at low temperatures. Thus, the 
density transform tends to be oscillatory and peaked at high wavenumber, so that large and finely-meshed grids are 
required. 

Fig. ^ shows the temperature of the marginal modes, presented in terms of the velocity dispersion and 
stability parameter Qs- The N — 2 results are probably not trustworthy to better than 2 s.f., but are included since 
they give at least an indication of the very low temperatures required. For = 3 and A^ = 4, the temperature 
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Figure 17. The dependence of the largest eigenvalue on growth rate and pattern speed for m = 3 disturbances in cut-out disks with 
N = 2 and /3 = 0.25. Note that the curves cross the real axis well to the left of the point (1,0) required for self-consistent modes. [The 
left-hand plot is for Qs = 1.0, i.e. cr^ = 0.283 and 7 = 11.0. The right-hand plot is for Qs = 0.5, i.e. 5-„ = 0.141, 7 = 48.5. The solid lines 
show curves of constant s at intervals of 0.1 from 8 = 0'*' (right-most curve) to s = 0.4 (left-most curve). The dotted lines show curves 
of constant Hp at intervals of 0.1 from Hp = 0.5 (lowest curve) to f2p = 1.3 (highest curve)]. 
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Figure 18. Minimum temperature for stability against m = 3 modes plotted against rotation curve index /3 for disks with various 
cut-out functions. The left-hand plot shows the minimum velocity dispersion 5-„; the right-hand one the same data presented in terms 
of the stability parameter Qs- The solid lines indicate the results for inner cut-out disks with Af = 2, 3, 4. 



is higher and hence the convergence is better. The rotation curve and inner cut-out affect the stabihty to three- 
armed perturbations in the same way as they affect the m = 2 stabihty. Disks where the centre has been cut out 
more sharply are more unstable than those where it has been removed relatively gently. Power-law disks with rising 
rotation curves are less securely stable than those with falling rotation curves. 



5.2.2 Neutral Modes 

So far, it appears as if the stability of the power-law disks to m = 3 perturbations is much like that to m = 2 
perturbations, the only difference being that even lower temperatures are required before growing modes can be 
excited. This conclusion seems to be reasonable enough in the light of eigenvalue curves such as Fig. ^ However, a 
different picture emerges when we pursue the mathematical eigenvalues down to vanishingly low growth rates and 



pattern speeds. Figs. 19 and ^ show the dominant eigenvalue curves as s and tend to zero, for Qs = 1 and 
/3 = ±0.25. For f3 = 0.25, the behaviour is similar to that observed for m = 2 perturbations. As the pattern speed 
is decreased, the eigenvalue curves continue to arc clockwise around the origin, tending to a a constant value (for 
given growth rate) in the limit fip — s- 0. They appear to be tending towards the negative real axis in the twin limits 
s ^ 0, ilp 0. However, for (3 — —0.25 the behaviour is strikingly different. Fig. ^ is much closer to the m — 1 
eigenvalue curves shown in fig. ^ It indicates that that modes with very low pattern speeds and growth rates are 
possible even at the high temperature of Qs — 1.0. It appears that there is a dramatic change in the behaviour of 
the eigenvalue curves between /3 — -f 0.25 and /3 — —0.25. 

It has proved hard to investigate the precise value of /3 at which such neutral modes become possible, since for 
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Figure 19. The dependence of the largest mathematical eigenvalue on growth rate and pattern speed for m = 3 disturbances in a disk 
with N = 2, Qb = 1 and fS = 0.25. The horizontal axis is Re[A] and the vertical axis is Im[A]. [The solid lines show curves of constant s 
at 11 values of s: s = 10"'^, 0.0001, 0.0005, 0.001, 0.005, 0.01, 0.05, 0.1, 0.2, 0.3, 0.4. The dotted lines show curves of constant Up at 20 
values of Qp-. Qp = 0.0001, 0.0005, 0.001, 0.005, 0.01, 0.02, 0.05, 0.1, 0.2, 0.3 .. . 1.2, 1.3]. 




Figure 20. The dependence of the largest mathematical eigenvalue on growth rate and pattern speed for m = 3 disturbances in a disk 
with A'^ = 2, Qs = 1 and /3 = —0.25. The horizontal axis is Re[A] and the vertical axis is Im[A]. [The solid lines and dotted lines are as 
in Fig. Fgl. 



values of /3 slightly greater than —0.25, the convergence becomes problematic. Fig. ^ presents evidence which is at 
least consistent with the idea that P = —0.25 is the delimiting value. The left-hand plot of Fig. |l] shows marginal 
eigenvalue curves for seven positive values of f3. None admit neutral modes. The right-hand plot shows marginal 
eigenvalue curves for six values of /3 < —0.25. All admit neutral modes. Fig. |2l| is for disks at Qs — 1. Similar curves 
were obtained at Qs = 0.3 (not shown). The evidence suggests a link between the absence of neutral modes and 
the presence of orbits which close in an inertial frame. For (3 > —0.25, it appears that the presence of three-lobed 
closed orbits prevents the formation of these neutral modes, even at low temperatures. For /? < —0.25, neutral modes 
persist even at high temperatures. 



5.3 Tetraskele or Four-Armed Modes 

Let us begin our study by examining the eigenvalue curves at moderate pattern speeds. Fig. ^ shows sample curves 
for (3 — 0.25. As before, the left-hand plot is for Qs — 1, and the right-hand plot for Qs = 0.5. In the latter case, 
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Figure 21. Marginal eigenvalue curves for m = 3 disturbances in inner cut-out disks with N = 2, Qb = 1.0 and various values of 
13. The left panel shows these marginal eigenvalue curves for positive values of /3 > 0.2, namely /3 = 0.2, 0.25, 0.3 . . . 0.45, 0.5. The 
right panel shows the curves for values of [3 < —0.25, namely [3 = —0.25, -0.3, -0.35 . . . -0.45, -0.5. For intermediate values of (3, the 
marginal eigenvalue curves are tangled, and are not shown here. Their behaviour is closer to the curves with /3 > 0.2 than to those with 
f3 < —0.25. [Each plot shows the largest mathematical eigenvalues for vanishing growth rate, s = lO"'^, and 20 values of the pattern 
speed: Op = 0.0001, 0.0005, 0.001, 0.005, 0.01, 0.02, 0.05, 0.1, 0.2, 0.3, . . . 1.2, 1.3.] 




Figure 22. The dependence of the largest eigenvalue on growth rate and pattern speed for m = 4 disturbances in Af = 2 and (3 = 0.25 
disks. The left panel is for Qs = 1.0, the right panel for Qs = 0.5. [The solid lines show curves of constant s at intervals of 0.1 from 
5 = 0"*" (right-most curve) to s = 0.4 (left-most curve). The dotted lines show curves of constant Qp at intervals of 0.1 from Qp = 0.3 
(lowest curve) to Qp = 1.3 (highest curve)]. 



the eigenvalue curves are ragged, as different eigenvalues become dominant. It is immediately apparent that the 
power-law disks are very stable to to = 4 perturbations. Fig. shows how low the temperature must be, in terms 
of the velocity dispersion (t„ and the stability parameter Qs, before the marginal eigenvalue curve intersects the real 
axis at unity. Even lower temperatures are required than for m = 3 (c.f., Fig. |8|). At these low temperatures the 
convergence is delicate. The values shown for i7u,min and Qg in Fig. ^ are probably accurate to only 2 s.f. Once 
again, disks where the centre has been cut out more sharply are more unstable than those where it has been removed 
relatively gently. Disks with rising rotation curves are more unstable than those with falling rotation curves, both 
in terms of (TM,min, the absolute amount of velocity dispersion required, and in terms of Qs, the ratio of the velocity 
dispersion needed relative to that necessary for axisymmetric stability. 

In the absence of closed four-lobed orbits, neutral m — A modes may be possible. However, these are more easily 
vanquished by temperature than the neutral m — 3 modes. In disks with /3 > 2/9, the four-lobed orbits closed 
in the inertial frame seemingly preclude the existence of neutral modes with fJp = even at low temperatures. 
However, disks with f3 only slightly less than the dividing possess very slowly growing modes that persist to quite 
high temperatures {Qs = 1). So, their marginal eigenvalue curves end up on the positive real axis. 
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Figure 23. Minimum temperature for stability against m = 4 modes plotted against the rotation curve index /3 for cut-out disks. The 
left-hand plot shows the minimum velocity dispersion (Ju ; the right-hand one the same data presented in terms of the stability parameter 
Qb- The solid lines indicate the results for inner cut-out disks with A'^ = 1, 2, 3, 4. 



6 THE SELF-CONSISTENT DISK 

Our attention now turns to the self-consistent power-law disk. Here, there is no cut-out function removing part 



of the density, and the potential experienced by stars in the disk is entirely due to their own gravity. Section 6A 



considers whether growing modes are possible in the self-consistent disk, while Section 6^ demonstrates the existence 
of neutral modes. 



6.1 Growing Modes 

The self-consistent disk is perfectly scale-free and so possesses no intrinsic length- or time-scales. This leads to some 
rather strange effects. As originally pointed out by Kalnajs in 1974 (see Zang's 1976 thesis), all dependence on 
growth rate and pattern speed can now be factored out of the integral equation. If a normal mode solution to the 
integral equation is found for some growth rate and pattern speed, it can be scaled to provide solutions at all growth 
rates and pattern speeds. This is proved for the family of self-consistent, power-law disks in Appendix ^ We shall 
make the same point here by dimensional analysis. 

Let us recall (see e.g., Section 2 of Paper I) that the equilibrium of the self-consistent disk contains three 
parameters - namely, the reference radius Rq, the reference density Sq and reference velocity w^. They are related 
by w| = 2TrGT,oRoL{(3). At first sight, it thus appears that we have two "degrees of freedom". For example, we 
could choose Sq and Rq with the consequence that v/s is then fixed. If this were the case, Rq would provide a 
length-scale and Ro/vp a time-scale. However, one of these "degrees of freedom" is illusory. The surface density is 
E = J^oR^^^ R^^^^ , so we cannot tell the difference between disks with the same Eq^o"'^'^- We are therefore free 
only to choose Eq-Rq^'^j which then determines v^R^. The disk is characterised by only two quantities, one with 
dimensions [M][L]''~^, and the other with [L]^+''[T]~^. It is apparent that these cannot be combined to yield a 



quantity with dimensions of length or of time (Lyndcn-BcU & Lemos 1993). Of course, carving out the centre of the 



disk breaks the self-similarity. The radius i? = i?o at which the inner cut-out is applied defines a length-scale. A 
time-scale is provided by the period of a circular orbit at this radius. The Fredholm integral equation then depends 
on the frequency uo. We are free to adjust Q so as to obtain non-trivial solutions, or modes. The self-similarity 
of the scale-free disk removes this freedom. The only adjustable parameter is the temperature. If the disk is hot, 
no modes exist. Intuitively, one expects that as the velocity dispersion is diminished below a critical value, modes 
become possible. The symmetry of the disk prevents it from selecting modes with a special frequency (unless it is 
zero). Thus, the self-consistent disk may possess isolated neutral modes. As regards growing modes, it may admit 
everything (a two-dimensional continuum with all possible frequencies) or nothing. 



The Fredholm integral equation (BIC) is hard to solve numerically for the self-consistent disk, since the kernel is 



highly singular. It immediately seems unlikely that growing, non-axisymmetric modes could exist in the self-consistent 



disk. All the quantities in the Fredholm integral equation ( BIO ) are complex, so there are effectively two equations: 
Re(A) — 1 and Im(A) = must be simultaneously satisfied. In the cut-out case, there are two free parameters: the 
frequency uj and temperature 7. In the self-consistent case, there is only one free parameter: the temperature 7. It 
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Figure 24. The left panel shows the response function TZ{a,5-u) for axisymmetric disturbances in the /3 = 0.25 disk. TZ is plotted as 
a function of wavenumber a for different values of the velocity dispersion o-„. The dashed line is the curve with an = o'u.mim that is, 
the critical velocity dispersion for stability according to local theory. The right panel shows the marginal stability curve TZ = 1 for the 
/3 = 0.25 disk. The solid (dotted) lines indicate the result of the global (local) theory. [The temperature and wavenumber at which 
instability sets in are marked. These are niin = 0.281 for global stability, ct^ = 0.283 for local stability; Ou = 4.38 for global 
stability, «u = 4.55 for local stability]. 
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Figure 25. Minimum temperature of the self-consistent power-law disk for global axisymmetric stability plotted against rotation curve 
index f3. In the left-hand plot, this is presented as the minimum velocity dispersion CTu, and in the right-hand plot as the stability 
parameter Qs necessary to stabilise the disk to axisymmetric perturbations. The dotted line shows the minimum temperature as given 
by local stability. 



seems unlikely that both the real and the imaginary parts of the integral equation could simultaneously be satisfied 
at any temperature 7. It seems even less likely that this can occur over a range of temperatures 7, as must happen if 
the growing modes affect all disks below a critical temperature. In the light of our analysis of the cut-out disks, there 
is one further piece of evidence that can be adduced. Non- axisymmetric modes in disks are principally confined to 
the region between the inner and outer Lindblad radii. As the pattern speed becomes vanishingly small, the inner 
Lindblad radius moves well beyond the inner cut-out radius. Under such circumstances, the response of the inner 
cut-out disk must approach that of the self-consistent disk (c.f., Zang 1976). If a continuum of modes exists in the 
self-consistent disk, its presence should be sensed by the cut-out disks as the pattern speed is brought to zero. We 
would expect the eigenvalues to draw together and become virtually independent of growth rate as the pattern speed 
approaches zero. There is no evidence that this happens. Here, the reader is urged to look back at graphs such as 
Fig. ^ (m = 1), or Figs. |l^ - (m = 3). Altogether, the weight of the evidence seems to us to suggest that there 
is no continuum and that the self-consistent disks admit no growing non-axisymmetric modes at all. 

The case of the axisymmetric modes is different. The work in Section 3.5 of Paper I shows that the kernel 
of the integral equation is Hermitian and the eigenvalues are real and positive for the m = modes. This means 
that the imaginary part of the Fredholm integral equation vanishes identically. Thus, there is just one equation 
and one free parameter, which can in principle be adjusted to obtain growing axisymmetric modes. As the disk 
cannot distinguish between modes with different growth rates, the self-consistent disks can admit a one-dimensional 
continuum of growing axisymmetric modes. 
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6.2 Neutral modes 

For the self-consistent disks, the neutral modes must have not just s = but ilp = as well. As noted by Lynden-Bell 
& Lemos (1993), the neutral modes are pure logarithmic spirals. In this case, the transfer function, which describes 
how the self-consistent disk responds to forcing, can be streamlined from the cumbersome expression given in Section 

3.3 of Paper I to a simple delta-function. Let us write 5m(a, a') = TZS{a — a'), where TZ defines the response function. 
We derive helpful expressions for the response function in Appendix ^ and apply these formulae in the following 
two sub-sections that study the axisymmetric and bisymmetric neutral modes. 



6.2.1 Axisymmetric Neutral Modes 

The only growing modes admitted by the self-consistent disk are axisymmetric ones. These set in through the neutral 



modes (Lynden-Bell & Ostriker 1967; Kalnajs 1971). The form of the response function (C9) in this instance is shown 
in the left panel of Fig. Now, breathing modes are long wavelength instabilities with no radial nodes. If present, 
the response function crosses the self-consistent line 7?, = 1 as a ^ 0. Fig. p3 immediately makes clear that the stellar 



power-law disks do not admit breathing modes - in contradistinction to the gas power-law disk ( Lemos et al. 1991 ) . 
Breathing modes are analogous to the instability of stars with 7 < 4/3, where 7 is the ratio of the heat capacity at 
constant pressure to the heat capacity at constant volume (see Binney & Tremaine 1987, chap. 5). Breathing modes 
occur in gas disks because gas molecules have internal degrees of freedom, which absorb energy released in changes 
such as contraction. This energy is then not available to the translational degrees of freedom which contribute to the 
pressure. There is thus insufficient pressure to resist gravitational collapse. Conversely, when a stellar disk contracts, 
all the energy released contributes to increased random velocity, which tends to counteract the effect of the increased 
gravity. Thus stellar disks are much more stable than gas disks to breathing modes. 

The marginal stability curve 7?, = 1 is plotted in the plane of dimensionless wavenumber a and velocity dispersion 
(T„ in the right panel of Fig. |2j. This compares the marginal stability curves obtained using local theory (dotted 
lines) and global theory (solid lines). The agreement is very good, especially close to the critical temperature. Fig. 
shows the minimum temperature needed for axisymmetric stability both in terms of the velocity dispersion and 
the stability parameter Q^. It is clear in this plot that the local and global results coincide most closely at /3 = 0, 
where the rotation curve is flat. This is probably because Toomre's local analysis assumes that the stellar disk has 
a Maxwellian velocity distribution. This is an excellent approximation for the Toomre-Zang disks, less good for 
the other power-law disks. Local theory becomes discrepant only when /3 < —0.25. Toomre's (1964) local analysis 
assumes that the velocity dispersion is small compared to the circular velocity, and that the wavelength of any 
disturbance is small compared to the radius R at the point under consideration. In the present notation, these 
conditions are it„ ^ 1 and a 3> 1. They become harder to justify when /3 < —0.25. As P becomes more negative, 
the wavelength of the unstable perturbation becomes larger, undermining local theory. The local calculation then 
underestimates the velocity dispersion needed for axisymmetric stability. Power-law disks with rising rotation curves 
remain globally unstable even when they are hot enough to be locally stable. 



22 A^. W. Evans and J. C. A. Read 




Figure 27. Velocity dispersion CTu and stability parameter Qs plotted against rotation curve index /3 for an m = 2 modes in self-consistent 
and cut-out disks. The curves labelled with the inner cut-out index TV show the (t„ and Qs at which the inner cut-out disks are 
marginally stable to bisymmetric disturbances. The curve labelled "scale-free" shows the (Tu,min and Qs at which the self-consistent disk 
just fails to admit any bisymmetric neutral modes. 



6.2.2 Other Neutral Modes 

Let us close with a brief examination of the bisymmetric neutral modes admitted by the self-consistent disk. The 



response function 7?. for these modes is given by eq. (C6). Fig. |26| shows the dependence of the response function on the 
wavenumber a for different temperatures. It is difficult to obtain an accurate value for TZ when both temperature and 
wavenumber are high. In Fig. the curve for each temperature has been truncated when the inaccuracy becomes 
severe enough to be visible on the graph. When the temperature is low enough, the disk admits neutral m = 2 
modes. As the disk is heated, the TZ-a curve drops, until a critical temperature is reached at which TZ = 1 only 
at some wavenumber. Above this temperature, the response remains less than unity for all wavenumbers and there 
are no neutral modes. Fig. ^ demonstrates that disks with falling rotation curves admit neutral modes at higher 
values of the stability criterion Qg than disks with rising curves. For f3 = 0.25, the curve with Qg = 0.4 intersects 
the line TZ — 1, whereas for /3 = —0.25, this curve never reaches TZ = 1. The dependence of the critical temperature 
for neutral modes on f3 is shown in Fig. ^ It has the opposite behaviour to marginally-stable modes in the cut-out 
disks, which are shown on the same plot. 



7 CONCLUSIONS 

The work in this paper concludes a linear normal mode analysis of a family of self-gravitating, differentially rotating 
stellar dynamical disks in which the circular velocity varies as a power of the Galactocentric radius, i.e., ticirc oc 
Our analysis includes both self-consistent disks, in which the surface density and the potential are related through 
Poisson's equation, and cut-out disks, in which stars close to the centre (and sometimes also at large radii) are carved 
out. The analytic and numerical procedures have exploited the self-similar force field of the models to develop a 
fast code that finds the pattern speeds and growth rates of the normal modes. Some tables of the fastest growing 
modes are provided in Appendix A comparison of these linear stability results with modes as detected in N-body 
simulations is a worthwhile and interesting task. 

(1) The inner cut-out is crucial in allowing growing non-axisymmetric modes in the power-law disks. Physically, 
it is helpful to picture the instabilities as standing wave patterns caused by the interference of travelling waves in 
a resonant cavity. For m > 1, the inner Lindblad resonance absorbs incident waves and damps any disturbance. 
Growing modes have pattern speeds which are large enough to push the inner Lindblad resonance within the shield 
of the inner cut-out. Waves cannot propagate easily within the core where there is little active mass. The inner cut- 
out acts as a reflective barrier. So, it has a destabilising effect if it shields waves from the absorbing inner Lindblad 
resonance. The inner cut-out is most effective as a reflector of waves when it is abrupt. Disks with lower cut-out 
indices allow more of the waves to seep through to the resonance and are thus more securely stable. The = 1 disk 
has the mildest cut-out and is the most stable of all. 

(2) The stability to different harmonics is summarised in Fig. |2^, which shows the minimum velocity dispersion 
needed for global stability. For an inner cut-out index A'^ < 3, the disks are more stable to bar-like or bisymmetric 
modes than to axisymmetric disturbances, whereas for A^ > 3, the converse is true. Modes with m — 3 and rn = 4 
are always very hard to excite. The disks must be extremely cold before they admit growing modes with more than 
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two arms. The most serious instabilities to which all the cut-out disks are susceptible are one-armed instabilities. At 
outset, let us concede that our analysis permits the m = 1 modes to shift the barycentre of the galaxy (c.f. Zang 1976). 
This artefact arises from the rigid core which has been carved out of the galaxy. A full treatment should allow the 
core to move in response to the growing mode in such a way as to cancel the shifting of the barycentre. Nonetheless, 
we do not believe this to be a serious flaw - first, because the barycentre moves only slightly within the linear regime 
and second, because N-body simulations by Earn (1993) , using the code described in Earn fc ScUwood (1995) , appear 
to corroborate the conclusion that one-armed instabilities are the most serious. The prevalence of the one-armed 
modes can be easily understood at a simple level, as there is no inner Lindblad resonance for to = 1. This removes 
a powerful stabilising influence. 



What is the physical mechanism generating these one-armed modes? They are not akin to the edge modes (Er- 



ickson 1974; Toomre 1981; Sellwood 198£). There is typically only one edge mode characterised by a single growth 
rate and pattern speed. By contrast, the cut-out disks exhibit an embarrassingly large number of different one- 
armed unstable modes. Second, the one-armed modes are not generated by refraction of short trailing waves into 
long trailing waves (the 'WASER') even though the Q profile of the cut-out disks does rise steeply in the inner 
regions in the manner envisaged by Lin and co-workers (Mark 1976; Lau et al. 1976). A third possibility, suggested 



to us by Sellwood (1997, private communication), is that the cut-out may be essentially invisible to long- wavelength 
disturbances, thus allowing trailing waves to leak through the centre and emerge as leading waves with devastating 
consequences. This does not seem to be the case, as the dominant wavelengths of the unstable modes are not large 
compared to the cut-out radius. Rather, the numerical evidence particularly of Fig. |l^ seems to us to suggest that 
the trailing waves are refiected as leading waves at the inner cut-out. It is this that completes the feedback for the 
swing- amplifier and causes one-armed mayhem in the disk. 

(3) For all azimuthal wavenumbers, the unstable modes persist to higher temperatures and grow more vigorously 
if the rotation curve of the power-law disk is rising (/3 < 0) rather than falling (/3 > 0). One way of understanding 
this result is to recall that, even in local axisymmetric theory, the power-law disks with rising rotation curves are 
more susceptible to Jeans instabilities and their ilk. In the absence of pressure, they remain vulnerable at longer 
lengthscales. They require larger velocity dispersions (relative to the local circular speed) for stabilisation at short 
scales. For example, on moving from /3 = 0.5 to (3 = —0.5, the typical densities increase by a factor of roughly three 
(c.f.. Table 1) and this in itself indicates that instabilities of whatever sort are much more worrisome in the models 



with rising r otation curves {(3 < 0). A number of authors (Lovelace & Hohlfield 1978; Papaloizou & Savonjie 1991 



CoUett 1995 ) have pointed out that gradients and especially turning points in the ratio of vorticity to surface density 
(or potential vorticity can provoke instabilities. In the power-law disks, the potential vorticity behaves like 



c = 



n 



(16) 



The gradient of the potential vorticity actually vanishes in the Toomre-Zang disk {(3 = 0). In the more nearly 
Keplerian disks (/3 > 0), swing amplification is boosted by this effect. In the more nearly uniformly-rotating disks 
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{(3 < 0), swing amplification gets weakened. Somewhat surprisingly, the effect of gradient terms in the potential 
vorticity actually runs contrary to the numerical stability results. Any effect, though present, must be quite weak. 
This has been demonstrated by Alar Toomre in unpublished work. 

(4) The criterion suggested by Ostriker and Peebles (1973) for stability against bar-like modes is untrustworthy. It 
gravely overestimates the energy in the form of random motions necessary to achieve stability against bisymmetric 
distortions. The power-law disks are much more stable than the criterion suggests. Similar qualms have already been 



reported by others (Aoki et al. 1979; Toomre 1981) 



(5) There are no length-scales and time-scales in the scale- free disks. If any mode is admitted at some pattern speed 
and growth rate, then it must be present at all pattern speeds and growth rates. Here, the analysis falls short of 
a complete proof, but our belief is that such a two-dimensional continuum of non-axisymmetric modes does not 
exist. The weight of the evidence suggests that there is no continuum and that the self-consistent disks admit no 
growing non-axisymmetric modes whatsoever. There are two telling pieces of evidence in favour of this conclusion. 
The first is the non-Hermitian nature of the kernel in the Fredholm integral equation for non-axisymmetric modes, 
which coerces us into satisfying (if we can) two simultaneous integral equations with only one free parameter at 
our disposal. Second, the behaviour of the mathematical eigenvalues of the cut-out disks in the limit of vanishing 
pattern speed must mimic that of the self-consistent disks and it does not suggest in the slightest the existence of 
a continuum modes. At first sight, this conclusion is astonishing - as it applies even to the completely cold disk! 
But, without reflecting boundaries such as provided by the cut-outs, there is no resonant cavity and no possibility 



of un stable normal modes. Such an outco me is already familiar in the theory of gaseous accretion disks ( Papaloizou 
& Sa fonjie 198^ |Papaloizou fc Lin 1995| ) 



The case of the axisymmetric modes is somewhat different. Now, the kernel is Hermitian and the eigenvalues 
must be real and positive. There is only one integral equation to be solved with our single disposable parameter. The 
temperature marking the onset of the neutral modes is accurately given by local theory. On physical grounds, these 
neutral modes are followed by growing modes as the temperature of the disk is lowered further. But the self-similar 
disk cannot distinguish between modes with different growth rates. Below the critical temperature the self-consistent 
disks must admit a one-dimensional continuum of growing axisymmetric modes. 

(6) The local theory of Toomre (1964) gives an excellent description of the global axisymmetric stability of the 
cut-out and the self-consistent power-law disks. The only possible axisymmetric modes are Jeans modes. Breathing 



modes do not occur in the stellar power- law disks, although they do occur in the gas disks ( Lemos et al. 1991 ). This 
is because gas molecules have internal degrees of freedom which can absorb energy released in contraction. This 
energy is not available to the translational degrees of freedom and does not contribute to the gas pressure resisting 
collapse. By contrast, when a stellar disk contracts, all the energy released is converted to random motions. 

(7) Both the self-consistent and the cut-out disks admit neutral modes. For the self-consistent disks, any neutral 
modes must have not just s = but also flp — 0. The existence of neutral modes seems to be related to the absence 
of m-lobed closed orbits in the frame rotating with pattern speed Qp. The cleanest example of this occurs for neutral 
triskele modes in the cut-out and self-consistent disks. For such disks with /3 < —0.25, there are neutral m — 3, 
fip = modes independent of temperature. These are the very disks for which three-lobed orbits closed in the 
inertial frame do not occur. This is straightforward to understand. The existence of any three-lobed closed orbits in 
the f2p = frame would cause any neutral pattern to evolve. 

NWE and JCAR are particularly indebted to David Earn and Jerry Sellwood for the loan of computer code that 
simulates the evolution of the cut-out power-law disks using the smooth-field-particle code described in David Earn's 
Ph. D. thesis. In particular, David Earn kindly re-ran some of the models in his thesis so that the dominant unstable 
modes could be checked against our linear stability code. We wish to acknowledge valuable conversations with James 
Binney, Jim CoUett, David Earn, Jeremy Goodman, Donald Lynden-Bell, Jerry Sellwood, Alar Toomre and Scott 
Tremaine - all of whose insights have benefitted us enormously. Both Jerry Sellwood and Alar Toomre offerred 
extensive and useful comments on earlier versions of the manuscript. Finally, NWE and JCAR thank the Royal 
Society and the Particle Physics and Astronomy Research Council respectively for financial support. 
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APPENDIX A: AXISYMMETRIC STABILITY 



This appendix briefly discusses the global stability of the cut-out power-law disks to axisymmetric modes. In the 
axisymmetric case, the marginal modes demarking the boundary between stability and instability are those of zero 
growth rate and pattern speed ( Lynden-Bell fc Ostriker 1967| ; Kalnajs 1971 ). In fact, the limit of vanishing growth 
rate gives rise to some numerical delicacies in the calculation of the kernel of the Fredholm integral equation. The 
expression for the cut-out angular momentum function Fim is derived in Appendix C of Paper I. For m = 0, it 
becomes 
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For very negative values of ?7, the e""^ term in (A.3) ensures that the limit is zero. But for positive 77, the expression 
for Fi„i fails to converge to any definite limit. This might be expected to present problems for the numerical 
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Figure Al. Dependence of the largest mathematical eigenvalue on growth rate s and temperature Qs for axisymmetric disturbances 
in inner cut-out disks with N = 2. [Curves are drawn for six different values of Q^. Qs = 0.6, 0.8, 1.0, 1.2, 1.4, 1.6. The Qs = 1 curve 
is drawn with a bolder line. For the /3 = 0.25 disk, Qs = 1 corresponds to a velocity dispersion of = 0.283, and for /3 = —0.25, to 
au = 0.509]. 
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Figure A2. Minimum temperature for global axisymmetric stability plotted against rotation curve index (3 for disks with various cut- 
out functions. The left panel shows the minimum velocity dispersion au\ the right panel shows the same data in terms of the stability 
parameter Q^. The dotted lines show the local theory; other lines indicate results of global theory for an inner cut-out disks. [Results for 
A'^ = 1, 2, 3, 4 are plotted in both graphs, but coincide so closely as to be virtually indistinguishable.] 



implementation. Fortunately, it turns out that the mathematical eigenvalue converges even though Fim does not 
(see Read 1997 for ample numerical evidence). This is because - especially when the growth rate is small - the 
eigenvalue is dominated by the matrix elements along the diagonal 77 = 0, where the troublesome term is simply 
unity. In eq. ( A.1), we can write N{1 — e~^^/^) Ki Trfj close to rj = 0. The dependence on the inner cut-out function N 
is cancelled out in the matrix elements close to the diagonal, where the transfer function is largest. We thus expect 
the stability of the inner cut-out disk to be virtually independent of - as already argued on physical grounds. 

Fig. Al shows how the mathematical eigenvalue depends on growth rate and temperature for two sample cut-out 
disks. For sufficiently small Qs, a unit eigenvalue can always be found for large enough growth rate s. When the 
disk is cool, it is susceptible to fiercely growing modes. As the disk is heated, the growth rate of the unstable mode 
decreases. The temperature at which the disk is marginally stable is given by the Qs for which the s = 0+ eigenvalue 
is unity. Inspection of Fig. Al shows that this value of Qs is close to unity. This already suggests that the global 
axisymmetric stability of the cut-out disks is close to the local axisymmetric stability. Fig. A2 shows the velocity 
dispersion CTu^min and the stability parameter Qs of the marginally stable disks as a function of rotation curve index 
p. Four different values of N are plotted, but the curves are scarcely distinguishable. For comparison, the results of 
local theory for the self-consistent disk are shown in a dotted line. Local theory is an excellent guide to the global 
stability of the inner cut-out disks. This is in accordance with Toomre and Zang's find ings - for t he disk with a flat 
rotation curve, local and global results agree "to within a few tenths of one per cent" ( Zang 1976 ). The value of the 
inner cut-out index has little effect on the stability. An additional outer cut-out function tends to make the disk 
slightly more stable, in that a lower velocity dispersion is needed for axisymmetric stability. 
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APPENDIX B: THE SYMMETRY OF THE SELF-CONSISTENT DISK 



In this appendix, we demonstrate that if any modes with uj ^ exist in the self-consistent disk, then they must form 
a two-dimensional continuum. This proof is a straightforward extension of that due to Kalnajs in 1974 for the disk 
with a flat rotation curve (given in Zang's (1976) thesis). 

For the scale-free disk, the angular momentum function is defined as 



2 + 27-/37 



Ik + mf2 — CjLJ^ ^ 



(Bl) 



where the notation of Paper I is used. As shown in Appendix C of Paper I, the angular momentum function consists 
of a regular function and a Dirac delta-function. Here, let us write: 



where 
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and t) = 2rj/{2 + /3). (For m 7^ 0, these expressions hold for all /. For m = 0, there is an ambiguity at the / = 



harmonic. The correct choices are F^^ 
delta-function in 77: 
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0). When lj — the angular momentum function is simply a 
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When = 0, we do not pick up the "off-diagonal" terms describing the response at wavelengths different from 
that of the original perturbation. Let us note that F^^lif) has no unique limiting value as w — *■ 0. It does not 
tend to zero in this limit and the case lo — Q cannot be recovered by taking the limit u; — > 0. Physically, this is a 
consequence of the symmetry of the disk. A pure logarithmic spiral is self-similar and a neutral logarithmic spiral 
introduces no length- or time-scales through its growth rate and pattern speed. Thus the response must be purely 
self-similar, i.e., another logarithmic spiral. Once the logarithmic spiral is allowed to grow or rotate - however slowly 
- a length- and time-scale has been introduced. The response is no longer self-similar, but involves contributions at 
many wavelengths. 

Analogously to (B2), let us divide the transfer hmction into a regular part and a Dirac delta- function: 

S,n{a,a) = S^^\a,a) + 5{a — a)S'^\a,a). (B6) 

sin {a, a') and Sm\a,a') may be obtained from eq. (78) of Paper I, using and respectively in place of 
Fim ■ We see that fI^ {rf) admits the factorisation 
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where Fi^{rj) is independent of the growth rate and pattern speed. The transfer function depends on the growth 

(2) 

rate only through the angular momentum function. So, Sm {a,a') is independent of Co, while the dependence of 
slri^{a, a') on uj can be factorised out as 

5W(a,a') - e«("-" )i"-5W(a,a'), (B8) 



where a is defined by a = 2a/(2 + /3) and sln' [a, a') is independent of Cj. Let us define reduced transforms via: 



A{a) = e 

Then, the Fredholm integral equation becomes 
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(B9) 
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where none of the quantities depend on uj. This means that if a non-trivial solution A{a) exists at all, then self- 
consistent modes A{a) exist with any non-vanishing growth rate and pattern speed. 
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APPENDIX C: NEUTRAL MODES OF THE SELF-CONSISTENT DISK 

In this appendix, expressions are derived for the response of the self-consistent disk to elementary forcing by a pure 



logarithmic spiral. For the neutral modes, the angular momentum function is simply a delta- function (B5). Thus, 
the transfer function is also proportional to a delta-function. Let us write Sm{ot,Oi') = TZS{a — a'), where TZ is the 
response function: 
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The response density transform Ares (a) at any wavenumber a is simply proportional to the imposed density transform 
Ahnp{a). Physically, the ratio of the response density to the imposed perturbation is given by the response function 
TZ. U TZ > 1, the perturbation grows; if 7?. < 1, the perturbation dies away. Self-consistent neutral modes require 
TZ — 1. We note that TZ is real and independent of the radial position in the disk, so that the disk is either stable 
or unstable everywhere. Recalling the symmetry properties of the Kalnajs function and of the Fourier coefficients, 
it is evident that the response function is independent of the sign of the wavenumber. This is an expression of the 
anti-spiral theorem (Lyndcn-BcU & Ostriker 1967). Every spiral of wavenumber a comes with an "anti-spiral" of 
wavenumber —a. Hence any neutral modes come in pairs, one leading and one trailing. 

There is an alternative form of the response function (CI), which merely involves integration over one period of 
the unperturbed orbit. This is especially useful for the important case of the m — neutral modes, as it resolves the 
ambiguity at the I = harmonic. In this derivation, the response to a pure logarithmic spiral is itself assumed to be 
a pure logarithmic spiral. According to Zang (1976), this calculation was carried out for the axisymmetric neutral 
modes of the disk with a completely flat rotation curve by Toomre in the 1960s. Putting together eqs. (51), (58) and 
(68) of Paper I, we obtain for the change in the distribution hmction 
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Substituting the potential of a single logarithmic spiral, we obtain 
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The integral in this equation can be split up into an infinite sequence of integrals over the radial period T . Summing 
the geometric progression, we have 
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For m 7^ 0, we simply set cj = to obtain 
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This expression can then be integrated over velocity space to obtain the response density Ercs- Recalling the definition 
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of the response function TZ, we finally obtain 
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This expression can be shown to be equivalent to that obtained directly from the Fredholm integral equation (see 
Read 1997). 

When the imposed perturbation is axisymmctric, cq. (04) is singular, viz.. 
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Using I'Hopital's rule to take the limit as s ^ 0, we find 
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This can be integrated over velocity to obtain the axisymmetric response function 
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In practice, these forms of the response function for neutral modes are more useful than (|C1| ) because they involve 
a quadrature performed over just one radial oscillation of the orbit. 



APPENDIX D: FASTEST GROWING MODES OF THE CUT-OUT DISKS 

The following sets of tables compare the stability of inner cut-out disks to perturbations of different azimuthal 
symmetry m. For each symmetry m and inner cut-out index A^, the pattern speed and growth rate of the fastest- 
growing mode are recorded. The results for (3 = 0.25, P = 0.00 and P = —0.25 are shown in different tables. The 
numerical accuracy parameters are: n = 301, Aa = 0.1, /min = —30, /max = +40, ncL = ^, fa — 0.8, Cacc = 20, 
ba.cc = 2, e\ — 10^^. The data in the tables are also presented in graphical form in the Figs Dl and D2. The pattern 
speed of the fastest-growing mode depends almost linearly on N. The fastest growth rate increases quite steeply 
with N, reflecting the increased instability of the more sharply cut-out disks. 
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Figure Dl. The variation of the growth rate and pattern speed of the fastest-growing m = 1 modes plotted against cut-out index N . 
All disks have Qa = 1.0. The solid line marked with diamonds is /3 = 0.25; the dotted line with triangles is /3 = 0.00; the dashed line 
with squares is ,8 = —0.25. 
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Figure D2. The variation of the growth rate and the pattern speed of the fastest-growing m = 2 modes plotted against cut-out index 
N. All disks have Qs = 1-0. The solid line marked with diamonds is (3 = 0.25; the dotted line with triangles is (3 = 0.00; the dashed line 
with squares is /? = —0.25. 



